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PREFACE 


If  experimental  work  is  to  be  accomplished  involving 
the  use  of  a  nuclear  reactor,  accurate  knowledge  concerning 
the  density  of  neutrons  at  various  locations  within  the  re- 
actor is  a  necessity,  whether  the  reactor  be  critical  or 
subcr i t i  cal . 

The  neutron  density  distribution  in  a  given  reactor 
can  be  estimated  mathematically  by  the  use  of  several  tech- 
niques including  transport  theory,  multigroup  diffusion  the- 
ory, and  single  group  diffusion  theory.   Of  these  various 
mathematical  theories,  the  single  group  diffusion  is  the 
only  one  amenable  to  hand  calculations.   Unfortunately,  how- 
ever, diffusion  theory  is  also  the  least  accurate  of  the 
three  named,  although  it  does  provide  a  good  first  approxi- 
mation to  the  actual  neutron  density  distribution  within  a 
reac  tor . 

In  order  to  accurately  determine  the  neutron  den- 
sity distribution  within  a  reactor,  experimental  measure- 
ments must  be  conducted  using  one  of  the  various  methods 
available  for  such  measurements.   One  of  the  most  convenient 
methods  is  that  of  foil  activation  measurements. 


Control  rods  made  of  materials  which  can  absorb 
great  numbers  of  neutrons  are  used  to  control  the  neutron 
distribution  within  reactors.   In  critical  systems.,  such 
control  rods  are;  of  course,,  necessary  for  safe  and  proper 
operation  of  the  reactor. 

The  University  of  Texas  maintains  a  subcritical 
reactor  which  is  used  for  teaching  purposes  and  for  research, 
The  purpose  of  this  thesis  is  to  accurately  determine  the 
effect  of  a  pure  cadmium  control  rod  on  the  neutron  density 
distribution  in  this  reactor. 

The  neutron  density  distribution  is  determined  by 
foil  activation  measurements  both  with  and  without  the  cad- 
mium control  rod.   Single  group  diffusion  theory  is  used  to 
predict  the  distribution  mathematically,  and  the  theoretical 
distributions  are  compared  with  the  experimentally  determined 
ones  . 

The  experimentally  determined  neutron  density  dis- 
tribution without  the  control  rod  is  compared  with  the  dis- 
tribution with  the  cadmium  control  rod  in  place  in  order  to 
clearly  show  the  effect  of  the  rod  on  the  neutron  density. 

The  author  wishes  to  express  his  appreciation  to 
Dr.  R.  N.  Little  for  his  advice  and  guidance  during  his 
supervision  of  this  research  and  to  Dr.  J.  A.  Scanlan  for 


his  many  helpful  suggestions  concerning  the  writing  of  the 
thesis.   Sincere  gratitude  is  also  due  Mr.  S.  C.  Mathur  and 
Mr.  C.  J.  Spahn  for  their  patient,  expert,  and  time  con- 
summing  assistance  with  the  equipment  used  in  this  experiment, 

Since  the  writing  of  this  thesis  marks  the  end  of 
a  long,  hard,  three  year  struggle,  the  author  feels  that  it 
should  at  least  have  a  dedication.   Therefore  this  writing 
is  dedicated  to  the  person  who  has  provided  the  encouragement 
and  the  understanding  which  is  so  necessary  in  this  type  of 

work --to  my  wife,  Patti. 
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CHAPTER 


THEORY 


A .   One  Velocity  Group  Diffusion  Theory  Approximation 
"without  Control  Rods 


In  the  single  velocity  group  diffusion  theory  it  is 
assumed  that  all  neutrons  in  the  system  being  studied  have 
the  same  velocity,  which  in  the  present  case  is  2200  meters 
per  second,  the  velocity  of  thermal  neutrons.   If  the 
neutron  population  is  described  in  terms  of  the  neutron 
density,  n(r,t),  where 

n(r,t)  =  number  of  neutrons  of  speed  v  per  unit 
volume  around  space  point  r  at  time  t; 

then  the  neutron  flux,  0(r,t),  is  expressed  in  terms  of 
n(r,t)  by  the  usual  definition 

0(r,t)  =  vn(r,t) 

where,  in  this  case,  all  values  are  for  thermal  neutrons 
and  v  is  constant  and  uniform  throughout  the  reactor. 

A  verbal  statement  of  the  neutron  balance  equation 
for  a  differential  volume  element,  dV,  includes  a  term  to 


R,  V.  Meghreblian  and  D.  K.  Holmes,  Reactor  Analysis 
(New  York:   McGraw-Hill  Book  Co0,  I960),  Chap.  5. 


account  for  the  net  gain  or  loss  through  the  boundaries  of 
the  volume  element,  and  is 


rRate  of  increase 
of  the  neutron 
density  in  dV 
about  r 


number  of 
neutrons  pro- 
duced in  dV 
per  unit  time-1 


.number  of 
neutrons  ab- 
sorbed in  dV 
■per  unit  timeJ 


rnet  number  of 
neutrons  escap- 
ing from  dV 

uper  unit  time 


This  statement,  in  mathematical  form,  becomes  the  equation 


|fdV  =  S(r,t)dV  -  Za0(r,t)dV  -  &(r,t)dV 


where  S(r,t)  is  called  the  source  function  and  denotes  the 
number  of  neutrons  produced  per  unit  time  in  dV,  and  the 
term  &(r,t)  is  called  the  leakage  function  and  denotes  the 
let  loss  of  neutrons  from  dV  due  to  diffusion  through  the 
boundari  es . 

In  diffusion  theory  the  leakage  term  becomes 

X(r,t)  =  -  DV20  (r,t) 


and  the  neutron  balance  equation  becomes,  under  steady  state 
conditions  for  which  n(r,t)  is  constant, 


DV20(r)  +  Z  0(r)  =  S(r) 

a 


In  this  equation  D  is  called  the  diffusion  coefficient  and 
is  defined  to  be 


X 


D  = 


tr 


3[Z.   ,  n  -  Z  (cos  0)    ] 
total     s        avg J 


where  X.^r  is  the  transport  mean  free  path  for  the  medium 
containing  the  volume  element   dV  for  which  the  equation  was 
derived.   Za  is  the  macroscopic  absorption  cross  section  for 
this  medium,  and  S(r)  is  the  source  term.   As  indicated  above 


tr 


tZtotal  ~  Zs(cos  0)avg] 


where  Z   is  the  macroscopic  scattering  cross  section  for  the 

medium,  (cos  0)avo-  =  2/3A  and  is  the  average  scattering  angle 

for  neutrons  involved  in  scattering  collisions  with  nuclei  of 

the  medium.   Z^    ^    .,  is  the  sum  of  Z   and  Z  • 
total  s       a 

The  net  neutron  current  through  a  given  area  lo- 
cated at  r  is  given  by  the  vector  equation 


J(r)  =  -DV0(r) 


The  partial  neutron  currents  through  an  area  are  defined  by 
the  equations 


j     (r)    =   i   0(r)    -    R    |V0(r)|cos    9 
+  42  ] 


and 


j_(r)    =    I   0(r)    +    R    |V0(r)|cos    9n 


where  0   is  the  angle  between  the  unit  vector  normal  to  the 

n  ° 

area  and  the  unit  vector  parallel  to  the  vector  V0(r).   The 

-> 
difference  of  the  partial  currents  is  equal  to   J(r).   In 

other  words , 


J(r)  =  j(r)  -  j  (r) 
+        — 


The  boundary  conditions  which  must  be  used  in  obtaining  a 
solution  for  the  above  diffusion  theory  neutron  balance  equa- 
tion are  as  follows: 

(1)  $(r)  must  be  finite  and  nonnegative  in  all  regions 
of  the  medium  in  which  the  diffusion  equation  ap- 
plies . 

(2)  At  any  interface  between  two  media  with  different 
diffusion  coefficients,  the  partial  currents  must 
be  continuous  through  the  interface.   If  there  is 


no  external  source  at  an  interface,  then  the  neutron 
flux,  0(r),  and  the  net  current,  J(r),  must  be  con- 
tinuous at  the  interface. 
(3)  The  flux  0(r)  must  vanish  at  some  finite  distance 
outside  the  boundary  or  else  the  partial  current 
from  the  area  outside  the  boundary  must  be  zero. 
In  diffusion  theory,  the  distance  beyond  the  re- 
actor boundary  at  which  0(r)  =  0  is  called  the  ex- 
trapolation distance,  d  (=  0.7104  A.  ,  )  . 

Two  assumptions  were  made  in  this  treatment  con- 
cerning the  reactor  assembly  without  the  control  rod.   It 
was  assumed  that  the  air  gap  along  the  central  axis  of  the 
core  did  not  exist.   This  air  gap  actually  represented  only 
approximately  1.20  per  cent  of  the  volume  of  the  core;  so 
this  assumption  was  not  considered  unreasonable.   It  was 
also  assumed  that  a  point  source  of  neutrons  was  located  at 
the  center  of  the  reactor.   This  point  source  is  treated  as 
a  boundary  condition  in  the  following  mathematical  develop- 
ment . 

In  the  case  of  finite  cylindrical  geometry  the 
diffusion  equation  can  be  written 


-  DV20(p,z)  +  Ea0(p,z)  =  3(p,z) 


6- 
where  S(p,z)  =  vZr^>{p,z)    for  a  multiplying  medium  and  is 
zero  for  a  nonmul tiplying  medium.   Z„  is  the  macroscopic 
fission  cross  section  for  the  medium  and  v  is  the  average 
number  of  neutrons  resulting  from  a  single  fission  event. 
If  both  a  core  and  a  reflector  are  involved,  a 
diffusion  equation  must  be  written  for  both  regions  as  follows 

D  V2^c(p,z)  +  vZf  ^.(p,z)  -  Za  6{p,z)    =    0 


DRV2jMp,z)  -  Za  )Mp,z)  =  0 


a  J°R' 
R 


The  solutions  to  these  simultaneous  equations  must  satisfy 
all  the  boundary  conditions  previously  stated. 

The  above  equations  can  also  be  written 

V2^c(p,z)  +  B§^p,z)  =  0 


and 


V2^R(p,z)  -  kUA9,z)    =    0 


R^R 


where  B2  =  [vL    -  Z  l/D„  and  k2,  =  Z  /V-q-        In  neutron  dif- 
fusion  theory  a  term  called  diffusion  length,  L,  is  defined 
such  that 


L«  4 


a 


Using  this  definition  in  the  terms  above  gives 


2      1 
K.  .  =  

R    LS 

I 


and 


B 


2     _       1 
C      T2 


-  1 


c  L   a, 


k   -  1 

CO 


where  kM  =  vZ-p  /za   is  the  infinite  multiplication  constant 
for  a  medium  composed  of  the  same  material  as  the  reactor 
core  . 

In  cylindrical  coordinate  the  V2  operator  is 


V2  =  I  L_ 
P  dp 


A_ 


dz 


2  > 


so  the  above  diffusion  equations  become 


and 


I  d_ 

P  dp 

I  i_ 
p  dp 


I-  tAt>,z) 


b-  rc 


dz 


-  ^c(p,z)  +  B2^c(p,z)  =  0 


I-  jUp**) 


>P 


R 


+  ±-z     ^R(P,0   "   K^R(P,Z)   =   0 


both  of  which  can  be  solved  by  the  method  of  separation  of 
variables . 

For  the  core  it  is  assumed  that 


(p,z)  =  Rc(p)Zc(z) 


and  the  resultant  general  solutions  to  the  separated  dif- 
fusion equation  for  the  core  are,  if  B2  =  k2  +  k2, 

^  '      c     p     z ' 


Z  (z)  =  A   sink  z  +C   cos  k  z 
c   '     z       z      z       z 


r  (p)  =  A  J  (k  p)  +  C  Y  (k  p) 
c  v  M/     p  o  v  pM/     p  o  x  pM/ 


where  J  (k  p)  and  Y  (k  p)  are  zero  order  Bessel  functions 
of  the  first  and  second  kinds,  respectively.   The  function 
YQ(k  p),  which  approaches  negative  infinity  as  p  approaches 
the  value  of  zero,  "becomes  more  negative  for  small  values 
of  (k  p)  than  the  function  J0(kQp)  "becomes  positive.   Thus 
since  the  boundary  conditions  stipulate  that  ^(p,z)  can 
never  be  negative,  C   must  be  taken  to  be  zero.   Since 
sin  kzz  becomes  negative  when  z  becomes  negative,  this  term 
would  violate  the  physical  requirement  that  j^(p,z)  be  sym- 
metric about  the  z  =  0  plane.   Thus,  Az  must  also  be  taken 
to  be  zero. 

Hence,  the  solution  for  the  core  region  is 


(p,z)  =  A  J  (kp)  cos  k 


c  "  o  v  p ( 


For  the  reflector  it  is  also  assumed  that 


•R(p,z)  =  Vp)ZR(z) 


and  the  resultant  general  solutions  to  the  separated  dif- 
fusion equation  for  the  reflector  are,  if  k2  =  k2  -  k   , 

'      R     P     z  ' 


Z  (z)  =  A„  sin  k  z  +  C   cos  kz 
R        z       z      z       z 


Rr(p)  =  a0ioupp)  +  cqko(k0p) 


where  Iq(k:  p)  and  KQ  (  k  p  )  are  zero  order  modified  Bessel 
functions  of  the  first  and  second  kinds,  respectively.   Here 
again,  the  condition  of  symmetry  of  the  neutron  flux  about 
z  =  0  requires  that  k%    be  taken  to  be  zero. 

Hence,  the  solution  for  the  reflector  region  is 


^r(p,z)  =  [ARI  (k  p)  +  CRK  (k  p)]  cos  k  z 


'R'o^p 


R1^  V 


If  pQ  and  zQ  are  defined  to  be  the  extrapolated 
radius  and  height,  then  the  flux  must  be  zero  at  the  points 
where  p  =  pQ  or  where  z  =  zQ.,   This  boundary  condition  re- 
quires that  jZ$R(p,z)  =  0  when  p  =  pQ  or  when  z  =  zQ.   So 


Kzzo  =  "Z2  sln 


ce  cos 


jt/2  =  0.   Thus  k       =  it/ 2 


z       ~o 

This  boundary  condition  also  requires  that 


10 


[ARIoUpp)    +    cRK0Upp)]    =    0       when    P    =    PQ 


Thus 


A„    .    -    C„    MjW 


■B 


'R 


Io(KpPo) 


The  equations  for  the  flux  in  the  reactor  thus  become,  in 
summary, 


jZ$  (p,z)  =  A  J  (k  p)  cos  k  z 


J^R(p,z)  =  Cj 


Ko^Kpp)  ~ 


lo(^pPo)  I°UPp) 


COS 


2z. 


where  zQ  and  pQ  are  the  extrapolated  height  and  radius,  re- 
spectively. 

At  the  interface  between  the  core  and  the  re- 
flector, the  boundary  conditions  require  that  the  flux  and 
the  net  current  be  continuous  across  the  interface.   Since 
the  current  is  a  vector  quantity,  this  condition  requires 
that  the  components  of  the  current  vectors  be  continuous. 
Thus,  at  the  interface  denoted  by  height  Zj_  and  radius  p, 

J^c(p;Zi)  =  ^R(p,zi) 


and 


Jc  (p, z1)  =  JR  (p, z1) 


11 


which  leads  to  the  equations 


AcJo(k0p)  cos  kzZi  =  CR 


and 


DckzVo(kDp)  sin  kzzi 


K0Upp)  " 


kq(kqPq)  t  (k    n) 


cos 


2*o   X 


iiT  DRCR 


V  f      >   K°(ltpPQ)   T   /      >" 


o  v '  pMo 


s  in 


2  z    i 
^zo 


Dividing  the  second  of  these  equations  by  the  first  gives 
the  equation 


Kz    tan  Kzz±    -    2z^  D^  tan  2^  z±     , 


which  can  be  solved  by  numerical  methods  to  determine  kz 
after  the  values  of  D-^,  D  ,  z^  ,    and  zQ  are  computed  or  meas- 
ured for  the  reactor  being  studied. 

Once  the  values  for  k„  and  kz  are  known,  the  cor- 
responding values  of  k  and  k  can  be  computed  from  the  re- 
lations 


B?  =  K    +  kz 


and 


2       2       2 

KR  -  Ko    Kz 


12 

since  B2  and  k2  depend  on  the  physical  constants  of  the  re- 
c       R 

actor  and  these  constants  are  known. 

The  constant  A   is  related  to  the  constant  CD  by 
c  K 

the  interface  boundary  condition  that  the  flux  must  be  con- 
tinuous at  the  interface  denoted  by  radius  p.  and  the  height 
z.   Thus  <£>    (p.,z)  =  ^R(pn.  ,z)  or,  if  the  point  where  z    0  is 


RVMi 


taken, 


CR  =  A 


Jo<kO01> 


rK0(Kp0l)   -  .^°[KPP0)  Iq{k~) 


which  will  give  Cp  as  a  constant  multiplied  by  Ac. 

The  constant  A   can  be  evaluated  by  the  boundary 
condition  of  continuity  of  net  neutron  current  at  the  inter- 
face between  the  external  neutron  source  and  the  core.   If 
comparison  of  theory  with  experimentally  determined  flux 
distributions  is  desired,  then  the  maximum  value  of  the  flux 

can  be  used  for  A  „ 

c 

The  values  of  Dc,  Dp,  Bc,  k„,  kz,  Kp,  k    ,     kz,  and 
Cp  for  the  reactor  used  in  this  experiment  are  computed  in 
Appendix  A. 

If  the  neutron  density  distribution,  n(p,z), 
rather  than  the  flux  distribution,  ^(p,z),  is  desired,  the 
the  constants  A   and  Cp  are  simply  divided  by  the  thermal 


13 


neutron    velocity    (220,000    cm/sec)    to    obtain    the    equations 


n     (p,z)     =    A'J     (k    p)     cos    k    z 
c  copM/  z 


nR(p>z)      =     C 


V      (  \  K0(KPP0}       T      (  ^ 

K°UPP)     "    I     (k    p     )     IoUPP) 
O        p     o 


COS      KZZ 


where    A'     =    A_/vth    and    Cft    =    CR/vth. 


R  -  °R/ vth 


One  Velocity  Group  Diffusion  Theory 
Approximation  with  Control  Rods^ 


Since  the  diffusion  theory  differential  equations 
are  the  same  for  the  core  and  reflector  regions  regardless 
of  the  presence  of  a  control  rod,  the  general  solutions  for 
these  equations  remain  the  same  as  for  the  case  when  no  con- 
trol rod  is  present.   The  effect  of  the  control  rod  is  ac- 
counted for  by  the  boundary  conditions  used  to  evaluate  the 
constants . 

The  general  solutions  for  the  diffusion  equations 
are 


(p,z)  =  [AcJ0(kp)  +  0cYo(knp)]cos  k, 


c"o  v^p 


J^R(p,z)   =   tARIo(Knp)   +  CRK0(K-p) JCOS   K 


R^o^V 


R.  V.  Meghreblian  and  D.  K.  Holmes,  Reactor  Analy- 
sis (New  York:   McGraw-Hill  Book  Co.,  I960),  Chap.  11. 


14 

where    B2    E    k2    +    k2    and    k2    e    kp    -    k2    as    before,    and    J^Ck^p) 
cpz  Rpz  '  oNpM/ 

and    YQ(k    p)    are    zero    order    Bess  el    functions    of    the    first 


a 


nd  second  kinds,  respectively,  and  I0(kqp)  and  KQ ( k  p ) 
are  zero  order  modified  Bessel  functions  of  the  first  and 


second  kinds,  respectively. 

The  boundary  conditions  which  were  applicable  for 
the  case  with  no  control  rod  remain  applicable,  and  for  the 
case  of  a  black  control  rod  fully  inserted  along  the  axis 
of  a  cylindrical  reactor,  the  additional  boundary  condition 
that  ft c  =  0  at  the  point  where  p  =  bQ  must  be  imposed.  The 
radius  bQ  is  defined  by  the  relation 

bQ  =  b  -  e 

where  b  is  the  radius  of  the  control  rod  and  €  is  the  extra- 
polation distance  into  the  rod. 

When  a  black  control  rod  is  present,  the  flux  be- 
comes zero  for  p  >>  0 ,  and  so  the  term  YQ(k  p)  need  not  be 
eliminated  from  the  equation  for  ^  (p,z).   The  boundary 
condition  that  ^c(b0,z)  =  0  yields  the  relationship 


AcJo(Vo)  =  "  CcYo(Vo)- 


This  relation  can  be  used  to  eliminate  C   and  yield  the 


15 


equation 


c(p,z)  =  Ac 


J  (k  b  ) 

J°(kp0)  -^ivv 


cos  k  z  , 


where  Ac^k  ,  and  kz  are  yet  to  be  determined. 

The  boundary  condition  that  ^p(p,z)  -  0  at  the  ex- 
trapolated radius,  pQ,  and  height,  zQ,  yields,  as  for  the 
case  of  no  control  rod,  the  relationships 


z    2z 


and 


AR  =  -  C 


Ko( KpPp) 
R  I0UpPo> 


where  k   can  be  determined  from  the  definition  Kp  =  k^  -  Kf 
when  zQ  is  known. 

In  summary,  the  equations  for  the  flux  can  now  be 
written 


Mp,z)  =  A< 


J0 (k_bn ) 

J°(V>  -  Y  (A  )  Y°(kp"' 
o  v  p  o  y 


cos  kzz 


>R 


(p,z)  =  C 


R 


_  ,    v    K0(«pP0)  _  ,    s 

K°UPp)  "  I  (k  p  )  ^V0 
o  v  pMo ' 


COS   K7  Z 


where  kz,  bQ,  and  k   are  known. 
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The  "boundary  condition  for  continuity  of  flux  and 

net  neutron  current  are  used  to  determine  the  remaining  con- 

s  tants . 

At  the  interface  denoted  by  height  Zj_  and  radius  p, 

the  relations 


^(p.zi  )  =  *Mp;zi) 


and 


Jc   (p,  Z±)      =     JR   (pjZ-jJ 


must  be  true.   These  relations  again  lead  to  the  equation 


cos  k  z  . 
z  1 


cos  k    z . 
z  i 


Dckz  sin  kzZl    Drkz  sin  kzz± 


which  can  he  written 


kz  tan  kzz^  =  — —  — —  tan   rt   z^ 
Dc  2 z0      2 z0 


and  can  be  solved  by  numerical  methods  in  order  to  determine 

k  .   After  k   is  determined,  k   can  be  found  from  the  rela- 
z  z  >       p 

tion  B2  =  k2  +  k2,  which  has  previously  been  defined, 
c     p     z '  -^         J 

The  relationship  between  Ac  and  Cr  is  determined 
by  the  boundary  condition  that  <f)c    =  ^R  at  the  interface 
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denoted  "by  radius  p^  and  height  z.   This  condition  yields 
the  relationship 


°E  "  Ac 


J0U 0Pi)  - 


^i^ly  (X    P  ) 


Ko(kdd1)  " 


"k°(v°>~    ; 


if  the  z  =  0  plane  is  used  for  the  evaluation. 

In  a  similar  manner,  the  boundary  condition  that 

Jc   =  Jp   at  the  interface  denoted  by  radius  p^  and  the 

P      0 
z  =  0  plane  can  be  used  to  determine  Ac  since  this  would  give 

two  simultaneous  equations  containing  AQ  and  C^.   In  order 
to  correlate  theory  with  experiment,  Ac  can  alternatively  be 
taken  to  be  the  flux  observed  along  the  midplane  at  the  dis- 
tance p  for  which 


J°(V>  -  £{££7  Yo<V>> 

o   p  o ' 


=  1 


The  values  of  all  constants  except  Ac  and  C^  are 
the  same  as  when  no  control  rod  was  present.   The  values  of 
Ac  and  C^  for  the  reactor  used  in  this  experiment  are  com- 
puted in  Appendix.  A„ 
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If  the  neutron  density  distribution,  n(p,z),  is 
desired  rather  than  the  flux  distribution,  jz$(p,z),  the  con- 
stants Ac  and  C^  are  divided  "by  the  thermal  neutron  velocity 
(220,000  cm/sec)  to  obtain  the  equations 


n  (p,z)  =  A 


o  v  p  o  ' 


cos  kz z 


'R 


(p,z)  =  C£ 


K  U  p  ) 

K°(KPp)  ~  I  (K°n     )     Io(Kpp) 
ok  *ppo; 


COS   KZZ 


where  A£  =  Ac/vth  and  CR  =  CR/vth. 


CHAPTER  II 


NEUTRON    FLUX    MEASUREMENT    BY    THE    ACTIVATION    METHOD 


A,       General 

The  activation  method  provides  one  of  the  most  con- 
venient and  most  accurate  techniques  for  measuring  neutron 
flux.   This  method  is  "based  on  the  fact  that  many  elements 
become  radioactive  when  exposed  to  a  neutron  flux.   For  sev- 
eral elements,  this  activation  occurs  only  for  a  narrow 
range  of  neutron  energies.   By  measurement  of  the  induced 
activity,  the  neutron  flux  responsible  for  the  activation 
can  be  determined. 

To  measure  a  flux,  a  thin  foil  of  the  most  suitable 
material  is  first  irradiated  for  a  certain  time  interval  at 
the  desired  location.   The  foil  is  then  removed  from  the 
neutron  flux  and  the  induced  activity  is  measured  in  an  ap- 
propriate counting  device.   The  relationship  between  the 
measured  activity  and  the  neutron  flux  must  include  correc- 
tion terms  for  such  things  as  counter  efficiency  and  depres- 
sion of  neutron  flux  by  the  foil  being  used  to  perform  the 
measurement . 
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In  order  to  "be  suitable  for  this  type  of  measure- 
ment, the  foil  material  should  have  a  large  neutron  absorp- 
tion cross  section  for  the  energy  range  to  be  studied*   Also, 
the  decay  scheme  of  the  radioactive  isotope  produced  by  the 
neutron  activation  should  be  simple  enough  for  easy  measure- 
ment and  interpretation,  and  the  half  life  for  decay  of  this 
product  isotope  should  be  such  that  irradiation  and  count- 
ing periods  are  not  excessively  long. 

For  thermal  neutron  detection,  the  most  frequently 

3 

used  elements   which  have  strong  resonance  absorptions  and 
yield  induced  radioactivity  of  suitable  half  lives  and  radia- 
tions are  In,  Rh,  Ag.  Au ,  I,  and  Dy. 


B.   Use  of  Dysprosium  Foils  for  Thermal 
Neutron  Detection 


The  dimensions  of  the  subcritical  reactor  used  in 
this  project  made  necessary  the  use  of  a  Van  de  Graaff  ac- 
celerator to  produce  the  external  source  of  neutrons  neces- 
sary for  a  subcritical  reactor  (see  Chapter  III).,   Due  to 
the  operational  characteristics  of  the  accelerator  used,  it 
was  possible  to  maintain  a  constant  source  strength  for  only 


3C.  F.  Bonilla,  Nuclear  Engineering  (New  York 
McGraw-Hill  Book  Co.,  Inc.,  1957),  p.  129. 
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short  periods  of  time,,  and,  hence,  the  foil  detector  ma- 
terial used  had  to  he  such  that  maximum  induced  activity 
could  be  reached  in  the  shortest  period  of  time  feasible. 
Of  the  foil  materials  (in,  Au,  and  By)  available,  dysprosium, 
which  has  a  large  cross  section  for  thermal  neutron  absorp- 
tion as  well  as  a  low  cross  section  for  epithermal  neutron 
absorption,  best  met  this  requirement.   However,  dysprosium 
does  have  one  disadvantage  in  that  approximately  twenty 
hours  must  be  allowed  for  adequate  decay  of  the  induced  ac- 
tivity prior  to  reusing  a  particular  foil  for  another  meas- 
urement . 

Preliminary  calculations  showed  that,  of  the  sev- 
eral naturally  occurring  isotopes  of  dysprosium  which  could 
result  in  radioactive  products  when  irradiated  by  neutrons, 
only  the  isotope  Dy164  need  be  considered.   Contributions 
to  induced  activity  from  neutron  absorption  by  isotopes 
Dy15   and  Dy158  are  negligible  due  to  low  percentages  of 
abundance  and  small  neutron  absorption  cross  sections  as 
well  as  the  long  half  life  of  decay  of  the  radioactive  prod- 
ucts  Dy     and  Dyx   .   The  contribution  to  measurable  ac- 
tivity by  the  transformation  of  Dy165  into  the  radioactive 
isotope  Dy166  was  determined  by  calculation  to  be  less  than 
2  X  10~14  per  cent.   These  calculations  were  confirmed  by 
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experimental  determination  of  the  half  life  of  decay  for 
the  activity  induced  in  natural  dysprosium  by  neutron  ac- 
tivation.  The  observed  half  life  was  exactly  that  of  Dy165 

Thus,  only  the  activity  induced  through  neutron 
absorption  by  the  isotope  Dy164  need  be  considered  in  con- 
nection with  neutron  flux  measurements.   The  reactions  in- 
volved are: 


Dy 


164 


+  n  -*  Dy 


is  sm 


Dy 


165 


+  7 


Ho165  +  p   +7 


Dy164  +  n  -  Dy165  ->  Ho165  +  p"+  y 

The  decay  scheme  for  Dy     is  shown  in  Figure   II-l, 

The  general  differential  equation  for  formation 
of  radioactive  product  nuclei  by  the  neutron  absorption  re- 
action is,  in  word  form, 


Rate  of  accumula- 
tion of  radio- 
active nuclei 


Rate  of  produc' 
tion  of  radio- 
active nuclei 


Rate  of  loss 
of  radio- 
active nuclei 


.164 


In  the  case  of  Dy1    two  differential  equations  are  neces- 


sary since  Dy164  can  be  transformed  into  both  Dy 


16  5 


and 
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Dyi65m     (i.25m) 


0.108  Mev 

IT 

9  7.6$ 

(2.3h) 


(3  decay 
2..H 


Dy 


IS  5 


** 


(35    0.305 


$2    1.000 


PI    1.215 


30    1.305 


Mev 


0.  985 


0.555 
0.515 

0.356 
0.275 


0.  095 


16  5 


Figure  II-l 


Decay  Scheme  for  Dy165  (from 

Tornau*) 


*Rudiger  Tornau,  Zeitschrift  fur  Physik,  159 ,    p 


101,  I960 
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Dy165m  by  the  neutron  absorption  reaction.   These  equations 


are 


(1) 


A  Dr™  (T) 


5m 


dT 


^Vth  -  ^mD5ffl(T) 


(2) 


d  D5  (T) 
"dT 


a4~5D.^,,   +  0.976  \c  D^  (T) 

n    4^th  5m  5mv  ' 


-  svt>  -  %srD5(T^th 


The  notation  used  in  the  above  equations  is  as  follows 


D, 


=  number  of  Dy164  nuclei  present  (assumed  to  be 
cons  tant ) 


Dc(T)   =  number  of  Dy1    nuclei  present  at  time  T. 


D5m(T)  =  number  of  Dy1S5m  nuclei  present  at  time  T, 


T 


=  irradiation  time 


fif-u  =  thermal  neutron  flux  during  irradiation 


f4-5 


=  thermal  neutron  absorption  cross  section  for 
activation  of  Dy164  to  Dy165. 


a 


4-5m 
ny 


=  thermal  neutron  absorption  cross  section  for 


164 


activation  of  Dyxo*  to  Dy 


.16501 


r5"6 

ny 


=  thermal  neutron  absorption  cross  section  for 
activation  of  Dy165  to  Dy166. 
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=    decay    constant    for    Dy 


16  5 


5m 


decay  constant  for  Dy 


16  5m 


Preliminary  calculations  showed  that  the  effect  of 
the  activation  of  Dy165  to  Dy166  is  completely  negligible, 
so  the  solutions  to  these  two  simultaneous  differential  equa' 
tions  are  found,  "by  standard  methods,  to  he 


d«„(t>  -  -s* 


°t:*a»j 


4r  th 


5m 


\ 


5m 


1  -  ex 


P  (-X5mT) 


and 


where 


D5(T)  =  B 


1  -  exp  (-\gT)   +  C 


exp  (-\_  T)  -  exp  (-\CT) 
5m  5 


(j4"5D/1  +  0.97  6  ct4"5^ 


B  = 


ny   4 


ny       4 


K 


fith 


and 


C  = 


0.976  a4~5mD, 
ny 4 

^5m  "  H) 


th 


In  determining  these  solutions  the  boundary  condition  that 

both  D,_(T)  and  D_  (T)  were  zero  at  time  T  =  0  was  used. 
5   '       5m 

If  the  foil  is  removed  from  the  neutron  flux  at 
time  T,  the  number  of  DylS5  and  Dy165   nuclei  present  at  any 
time  t  thereafter  is  determined  by  the  solutions  to  the  dif- 
ferential equations 
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*DKB1(t) 


5m 


dt 


and 


dD5(t) 


dt 


"  X5mD5m("fc) 


=  0.976  XHmDfim(t)  -  XKDK(t) 


5^5 


In  this  case  the  boundary  conditions  to  be  used  are  that 
D5m(t)  =  D5m(T)  and  D5 ( t )  =  D5(T)  at  t  =  0,  where  t  is  time 
after  removal  from  the  neutron  flux. 

The  solutions  to  these  simultaneous  differential 
equations  are 


and 


D.Sm(t)  -  DSm(T) 


5m 


exP  (-^Kmt) 


»] 


D5(t)  =  ([      -     X^f5^^^    ("  X5mt)"  6XP  ("  X5t) 


+  D5(T) Texp  (-  X5t) 


Since  the  activity  at  any  time  t  is  simply  the 
product  of  the  decay  constant  and  the  number  of  nuclei  pres- 
ent at  that  time,  the  following  equations  give  the  activity 
as  a  function  of  the  time  after  removal  of  the  detector  foil 
from  the  neutron  flux: 


A5m(t)  =  X5mD5m(t) 


A5(t)  =  *.5D5(t) 
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Thus,  if  the  activity  is  counted  for  a  specific 
time  interval,  t   to  t  ,    the  number  of  radioactive  disinte- 
grations, N,  occurring  during  this  time  will  be 


A5m(t)  dt  +  I     A5(t)  dt 


\c    D,  (t)  dt  +  I     X._D_(t)  dt 
5m  5m   '        I      5  5 


where  N  =  total  number  of  disintegrations  occurring  between 
times  t,  and  t„ .   Evaluation  of  this  integral  yields 


N  = 


0.976  X 


5m 


5     5m 


exp  (-  ^5ti)  -  exp  (-  A-5t2) 


1-976    X5    ~    X5m 


(Vr  "  \Rm) 


Dp-   (T) 

5m  ' 


exp  (-A-,-.  t,  )  -  exp  (-  A,.-  t„) 

*  5m  1'      *  5m  2' 


This  equation  for  N  appears  to  be  extremely  cum- 
bersome and  unwieldy,  but  after  the  constants  are  evaluated 
and  the  time  intervals  to  be  used  for  irradiation  and  count- 
ing are  decided  upon,  the  above  equation  reduces  to  N  = 
constant  X  m  $ .,     where  m  is  the  weight  of  the  detector  foil 
in  grams . 
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Trial  runs  and  calculations  using  foils  of  various 
weights  indicated  that  an  irradiation  time  of  five  minutes 
( T  =  0  to  T  =  5 )  and  a  counting  interval  of  10  minutes 
(t   =  15  to  t   =  25)  would  give  very  accurate  results.   Evalu^ 
ation  of  the  various  constants  needed  for  the  above  equa- 
tions was  carried  out  and  the  results  are  tabulated  in 
Table  II-l.   Using  these  constants,  the  equation  for  N  re- 
duces to 


N  =  0.6402  m^th 


where  m  =  foil  mass  in  grams  and  fi^    =  neutron  flux  in 
neutrons  per  cm2-min.   After  converting  the  flux  into  neu- 
trons per  square  centimeter  per  second,  the  equation  for 
thermal  neutron  flux  can  be  written  as  follows: 


/neutrons \   N 

thV  cm2  secy   38.412  m(grams) 


C .   Determination  of  Thermal  Neutron  Flux 
from  Measured  Activity 


The  measured  activity  must  be  multiplied  or  di- 
vided by  several  correction  factors  in  order  to  account  for 
losses  occurring  in  the  counting  process,  according  to 
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TABLE   II-l 

CONSTANTS   TO  BE  USED  IN  EQUATION  FOR   NUMBER   OF 

DISINTEGRATIONS  DURING  COUNTING  INTERVAL 


Cross   Sections* 


a4-5   =  800  X  10~24cm2 
n^ 


a4-5m  _   2000  x  10-24cm2 
ny 


Decay  Constants 


\tr    =       0.695        =  0.004978  min"1 


139.2  min 


X5m  =  0.695 =  0.5544  min"1 

1.25  min 


Product  of  Decay- 
Constants  and 
Times 

"T  =  5  min" 
tj  =  15  min 
to  =  25  min_ 


\5T  =  0.02489 


\5t1  =  0.07467 
\5t2   =  0.12445 


k5mT      =      2.77  20 


Wl  =   8'3160 
^5mt2  =  13-8600 


Values  of 
Exponentials 


-X   T 
e  5     0.97541 


5  X  =  0.92805 


-K   T 
e  sm     0.062536 


e"  5m  1  =  Oo 000245 


Act 


5U2 


0.88298 


e"  5m  2  =  0.000001 


e  o  1  „  e  5  ^   =  0.04507 


e^5mtl  _  e-^5mt2  =  0. 000 244 


*Hughes,  Magurno,  and  Brussels  Neutron  Cross  Sections , 
Supplement  No.  1  to  BNL  525.  U.  S.  Government  Printing  Office,  Washing- 
ton 25 ,   D.  C.,  January  1960. 
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Greenfield,  Koontz,  Jarrett,  and  Taylor.4   When  the  activity 
being  counted  is  beta  activity,  losses  are  experienced  due 
to  self-absorption  and  scattering  by  the  foil  material  as 
well  as  the  usual  losses  due  to  counter  efficiency,  flux  de- 
pression, geometry  factor,  etc. 

If  N  is  the  number  of  disintegrations  occurring 
during  the  interval  in  which  the  neutron  detector  foil  is 
counted,  and  C,,   is  the  total  number  of  counts  during  this 
interval  which  are  due  to  activity  induced  by  thermal  neu- 
trons only,  then  the  relationship  between  N  and  C,,  is 


C.,  F    (d) 

N  _  th   sp 

1.011  F   F   F   FT  F       '    (~dT 
7   g   e   bs   s e&s  c 


where  the  factor  1.01.1  is  a  correction  factor  to  take  into  ac- 
count the  Maxwellian  distribution  of  thermal  neutrons  as  well 
as  the  dependence  on  neutron  energy  of  the  <r  values  for  the 
detector  foil  material.  The  F  factors  are  correction  factors 
for  gamma  photon  counts  from  competing  modes  of  decay  in  dys- 
prosium, detector  geometry,  counter  efficiency,  back  scatter- 
ing, self  absorption  and  self  scatter  in  the  foil,  and  flux 

depression  by  the  foil.   The  factors  F    and  F   „    are  depend- 
*  •>  Sp       sa&sc 

ent  on  the  foil  thickness,  d. 


Greenfield,  Koontz,  Jarrett,  and  Taylor,  Nu- 
cleoni  cs ,  15 ,  No.  3,  p.  57,  March,  1957. 
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The  total  counts  observed,  C!   ,  which  are  due  to 

7       th 

activity  induced  by  thermal  neutrons  is  determined  by  takin; 
one  measurement,  CC(~[,  in  a  certain  location  with  a  foil 
wrapped  in  a  cadmium  cover,  and  by  taking  a  second  measure- 
ment, C ,  ,  -J  ,  in  exactly  the  same  location  using  the  same 
foil  with  no  cadmium  cover.   The  difference  of  these  meas- 
urements is  the  desired  value,  C^,   The  mathematical  ex- 
pression for  C-th  is 


Cth  "  Ctot    Ccd 


when  the  foil  material  is  dysprosium.   For  materials  which 

have  large  resonances  in  their  absorption  cross  sections  for 

neutron  energies  greater  than  thermal,  the  cadmium  covered 

foil  measurement,  Cc(^,  must  be  multiplied  by  a  correction 

factor  F  ^  to  account  for  epithermal  neutrons  absorbed  by 

the  cadmium  cover.   However,  in  the  case  of  dysprosium,5 

there  are  no  such  large  resonances,  and  thus,  since  C  ,  is 

only  about  1  per  cent  of  C,   ,  -,  ,  this  correction  is  not 
J  *  total' 

necessary.   Preliminary  data  taken  using  dysprosium  foils 
showed  that  omission  of  this  correction  factor,  F  .,  which 


1951 


'C.  W.  Tittle,  Nucleonics  ,  8_,  No.  6,  p.  5,  June, 
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as  yet  has  not  been  experimentally  determined  for  the  case 
when  the  detector  foil  material  is  dysprosium,  would  intro- 
duce an  error  of  less  than  0.07  per  cent  in  the  data  taken 
for  this  thes  is . 

Since  the  beta  counting  apparatus  used  (see  Chap- 
ter III)  had  4jt  geometry,  the  correction  factor  for  geometry, 

F  ,  is  equal  to  unity.   The  scintillator  container  dimen- 
g 

sions  were  such  that  all  beta  particles  were  captured  prior 
to  their  reaching  the  walls  of  the  container.   Thus,  the  cor- 
rection factor  for  back  scattering,  F,  ,  is  also  equal  to 
unity.. 

All  gamma  rays  emitted  during  the  10-minute  count- 
ing interval  are  coincidental  with  the  beta  particles  which 
are  being  counted,  so  the  factor  F  ,  which  compensates  for 
gamma  photons  emitted  by  decay  modes  other  than  the  one  being 
measured,  is  also  equal  to  unity.   The  15-minute  waiting  period 
prior  to  the  counting  interval  was  purposely  chosen  so  that 
the  number  of  gamma  rays  from  the  isomeric  transition  of 
Dy165m  to  Dyi65  (half  life  of  1#25  min.)  would  be  completely 
negligible  in  that  less  than  0.025  per  cent  of  the  activity 
present  due  to  this  decay  mode  at  time  t  =  0  remains  at  time 
t]_  (=  15  minutes). 
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The  factor,  ^ST),    which  compensates  for  the  depres- 
sion of  the  thermal  neutron  flux  by  the  foils  being  used  to 
conduct  the  measurements,  can  be  calculated  by  the  equation6 


Fsp(d)  =  1  +  Q^i  aR      if  R  «  Xt. 

Hr 


which  was  developed  by  Bothe.   In  this  equation  R  =  foil 


radius  and  A. 


X, 


,   = 7 b  a  \ =  the  transport  mean  free 

tr    1  +  (cos  0) 


avg 


path  in  the  reactor.   Xs  is  the  scattering  mean  free  path 

in  the  reactor,  and  (cos  0)     =  2/3A  where  A  is  the  mass 

avg 

number  of  the  scattering  nucleus.   The  term  a    depends  on  the 
foil  thickness,  d,  and  is  determined  by  the  equation 


a 


(d)  =  1  -  [exp  (-  ud)]  (1  -  ud)  +  u2d2E.(-ud) 


where  d  is  thickness  of  the  detector  foil,  u  is  the  thermal 
neutron  absorption  coefficient  for  this  foil,  and  E^  (-|j.d) 
is  the  exponential  integral  function7  value  for  the  foil. 


1951. 


C.  W.  Tittle,  Nucleonics  ,  9_,  No.  1,  p.  60,  July, 


TE.  Jahnke,  and  F.  Emde,  Tables  of  Functions, 
Dover  Pub.  Co.,  New  York,  1943. 
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Calculations  using  this  equation  were  compared 

with  experimental  values  by  Bothe,8  and  agreement  was  excel- 
lent even  for  cases  where  the  foil  radius,  R,  was  only 
slightly  less  than  ^fr-        Thus  this  equation  was  used  to  cal- 
culate FSp(d)  for  the  foils  used  in  this  project.   The 
values  calculated  are  tabulated  in  Chapter  IV. 

The  factor  Fe  accounts  for  the  efficiency  of  the 
beta  detection  apparatus  and  is  less  than  unity  even  for  a 
4it  geometry  detector  since  the  contributions  by  low  energy 
beta  particles  may  be  lost  in  the  tube  noise  (see  Chapter 
III).   The  factor  F   p   (d),  the  correction  for  self  absorp- 
tion and  self  scattering  of  beta  particles  by  the  foil,  is 
dependent  on  foil  thickness  and  has  not  been  determined  for 
dysprosium  foils.   By  coincidence  measurements  utilizing  the 
fact  that  gamma  rays  are  emitted  simultaneously  with  beta 
particles  in  the  decay  of  Dy165,  the  combined  correction 
factor   F   F   »   (d)   can  be  determined.   The  equations  used 
in  this  determination  are: 


C   =  Kfef„N     For  the  gamma  ray  counter 

Co  =  FeFsa(^sCF   N    For  the  beta  particle  counter 


1951 


8C.  W.  Tittle,  Nucleonics  ,  9_,     No.  1,  p.  60,  July, 


and  Cp_7  =  Kf^f^F^F 


ex  g-1-  e1-  sa&s  c 

16  5 
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,N  for  the  coincidence  counter  where 
N    number  of  Dy-1-00  disintegrations  during  the  counting  inter- 
val and  the  C's  are  the  various  observed  counts  corrected 
for  background.   The  factor  [K  fg  f  ]  accounts  for  the  geom- 
etry and  efficiency  of  the  gamma  detector  and  for  gamma  rays 
absorbed  prior  to  reaching  the  gamma  ray  detector.    Elimina- 
tion of  [K  fe  f g ]  yields  the  equation,  since  Fg  is  known  to 
be  unity,  N  =  CpC  /Cg   .   Thus,  if  Co  for  the  4 it  counter  is 
then  measured, 


Cft(47t  counter) 

— E -  f 

N 


e  Fsa&sc  = 


Cp.  Co  (4  it  counter) 
C7Co, 


Since  F       is  a  function  of  foil  thickness,  d,  coincidence 
sa&sc  '   ' 

counting  measurements  were  conducted  (see  Chapter  III)  for 
each  foil,  and  the  resulting  values  for  [F   F   o    (d)]  are 
tabulated  in  Chapter  IV. 

Thus,  when  all  correction  factors  are  included, 
the  final  equation  for  absolute  thermal  neutron  flux  for  a 
given  number  of  counts  observed  during  the  counting  interval 
t-^  (=  15  minutes)  to  tg  (=  25  minutes)  is  given  by  the  equa- 
tion 


t 


[C.  .  ,  -  C  ,]  F   (d) 

L  total     cdJ   spv  ' 

th    [1.011  Fe  Fsa&sc(d)][38.412  m] 
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where  m  is  the  foil  weight  in  grains  and  jz^,,  is  in  neutrons 
per  square  centimeter  per  second. 

If  the  result  is  desired  in  terms  of  neutron 
density,  n  (neutrons/cm3),  we  must  divide  ^,,  "by  the  thermal 
neutron  velocity  of  220,000  centimeters  per  second  to  get 
the  equation 


n  = 


£Ctotal  "  Ccdl  Wd) 


220,000  [1.011  Fe  Fsa&sc(d) ] [38.412  m] 


where  n  is  in  neutrons  per  cubic  centimeter  and  m  is  the  foil 
weight  in  grams. 


CHAPTER     III 


EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


A .   Reactor  and  External  Neutron  Source 

The  University  of  Texas  subcritical  reactor9  is 
composed  of  a  polyethylene -clad,  homogeneous  core  of  20  per 
cent  enriched  UO2  dispersed  in  high  strength  Marlex-50 
polyethylene ,  and  three  separate  reflector  assemblies. 
The  core  is  a  cylinder  approximately  10  inches  in  diameter 
and  approximately  14  inches  long,  assembled  from  eight 
homogeneous  polyethylene  discs  impregnated  with  U02.   One 
radial  and  one  axial  hole  of  approximately  1  inch  diameter 
penetrate  the  core  for  the  purpose  of  insertion  of  foils, 
neutron  sources,  etc. 

The  reflector  assemblies  are  made  of  polyethylene 
and  graphite.   The  polyethylene  reflector  is  composed  of 
two  concentric  shells,  each  approximately  three  inches 
thick,  so  that  a  total  reflector  thickness  of  either  three 


9F .  D.  Eichenbaum,  Lockheed  Nuclear  Products 
Reports  NR  57  and  58,  March,  1959,  Nuclear  Products  Divi 
sion  of  Lockheed  Aircraft  Corp.,  Marietta,  Georgia. 
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or  six  inches  may  "be  used.   The  graphite  reflector  is  ap- 
proximately ten  inches  thick.   All  three  assemblies  are 
constructed  in  the  form  of  annular  cylinders  with  end  plugs 
to  fill  the  annulus  so  that  the  core  can  be  completely  re- 
flected . 

During  experimentation,  the  core  is  inserted  into 
one  of  the  reflector  shells  and  the  appropriate  plug  is 
inserted  into  the  end  gap,  thus  providing  a  uniform  thick- 
ness of  reflector  at  the  ends  and  around  the  core.   When 
the  core  is  properly  positioned,  its  axial  and  radial  ports 
are  aligned  with  similar  ports  provided  in  the  reflector. 
Thus  foils,  control  rods,  and  neutron  sources,  appropriate 
to  the  experiment  being  conducted,  can  be  positioned  in 
the  port  s . 

Small,  cylindrical  plugs,  some  composed  of  core 
material  and  some  of  reflector  material,  are  provided  in 
sufficient  quantity  for  filling,  as  appropriate  to  the 
experiment  being  conducted,  the  voids  created  in  the  re- 
actor by  the  radial  and  axial  access  ports. 

In  this  experiment  the  6  inch  polyethylene  re- 
flector was  used.   A  cross  sectional  view  of  the  core  with 
the  6  inch  polyethylene  reflector  in  place  is  shown  in 
Figure  III-l.   The  dimensions  in  this  figure  are  in 
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centimeters  and  are  accurate  dimensions  rather  than  the 

approximate  figures  used  in  the  above  description.   Dimen- 
sions of  the  individual  core  element  discs  as  well  as  the 
various  constants  applicable  to  the  core  and  the  reflector 
are  tabulated  in  Appendix  A. 

The  cadmium  control  rod  used  in  this  experiment 
was  a  hollow  cylindrical  tube  86.2  centimeters  long  com- 
posed of  commercially  pure  cadmium.   This  tube  had  an  out- 
side diameter  of  2.62  centimeters  and  a  tube  wall  thickness 
of  0.09  centimeters. 

The  boron  control  rod  was  a  hollow  cylindrical 
tube  80.4  centimeters  long  composed  of  polyethylene  4  per 
cent  impregnated  with  boron.   The  outside  diameter  of  the 
tube  was  2.71  centimeters  and  the  tube  wall  thickness  was 
0.13  centimeters. 

Because  the  University  of  Texas  Subcritical  Re- 
actor is  subcritical,  an  external  source  of  neutrons  must 
be  used  in  order  to  maintain  a  thermal  neutron  flux  within 
the  reactor  under  steady  state  conditions.   Several  manu- 
factured Pu-Be  sources  with  a  constant  neutron  yield  are 
available  for  use  with  the  reactor.   Unfortunately,  how- 
ever, the  dimensions  of  these  "fixed"  sources  prohibit 
their  use  inside  the  reactor  when  one  of  the  control  rods 
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is  positioned  in  the  axial  hole.   Thus,  it  was  necessary 
to  use  the  Van  de  Graaff  accelerator  to  create  the  neces- 
sary external  neutron  source. 

The  Van  de  Graaff  particle  accelerator  available 
for  use  in  connection  with  this  experiment  was  manufactured 
by  High  Voltage  Engineering  Corporation  and  is  capable  of 
producing  deuterons  accelerated  to  energies  of  up  to  1  Mev . 
A  beryllium  target  was  bombarded  by  1  Mev  deuterons  from 
the  accelerator  and  neutrons  were  produced  by  the 
Be9(d,n)B10  reaction.   The  location  of  this  external  neu- 
tron source  within  the  reactor  is  shown  in  Figure  III-l. 

Beryllium  was  chosen  as  the  target  material  be- 
cause the  stability  characteristics  of  the  accelerator 
required  operation  at  deuteron  energies  between  0.8  and 
1.0  Mev,  and,  in  this  deuteron  energy  range,  the  Be9(d,n)B10 
reaction  afforded  a  better  neutron  yield  than  does  any 
other  X(d,n)Y  reaction.   Also,  since  beryllium  targets 
were  available  in  metal  disc  form,  the  problems  associated 
with  target  cooling  and  handling  were  greatly  simplified 
by  use  of  beryllium  rather  than  tritium,  the  only  other 
target  material  immediately  available  for  use  in  connec- 
tion with  this  experiment.   Target  cooling  was  accomplished 
by  means  of  a  stream  of  air  blown  directly  on  the  target 
assembly . 
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The  energy  distribution  of  neutrons  from  the 
Be9(d,n)B10  reaction  for  1  Mev  deuterons  (as  shown  in 
Figure  III-2)  shows  a  preponderance  of  neutrons  having 
energies  of  approximately  3.9  Mev  and  4.7  Mev.   The  angu- 
lar distribution10  curves  for  the  neutrons  produced  by 
this  reaction  indicate  only  slight  directional  preferences, 
and  measurements  conducted  to  determine  the  strength  of 
the  source  used  in  this  experiment  showed  that,  if  an  iso- 
tropic neutron  yield  were  assumed,  the  maximum  deviation 
due  to  directional  preference  would  be  less  than  6.5  per 
cent.   Thus,  it  was  assumed  for  the  purpose  of  this  experi- 
ment that  the  Be9(d,n)B10  reaction  functioned  as  an  iso- 
tropic source  of  neutrons. 

For  this  experiment,  the  source  strength  had  to 
be  determined  as  a  function  of  the  electrical  current  cre- 
ated by  the  action  of  the  beam  of  accelerated  deuterons 
produced  by  the  accelerator  and  striking  the  beryllium 
target.   This  was  accomplished  by  measuring  with  a  "long 
counter"  the  number  of  neutrons  produced  at  various  beam 
current  settings.   These  measurements  were  taken  at  various 
angles  relative  to  the  deuteron  beam. 


10A.  I.  Shpetnyi,  Soviet  Physics  JETP,  5,  p.  357, 
October,  1957  . 
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N(E  )   (Summed  over 
all  angles  ) 


En  (Mev) 


Figure  III-2 


Energy  Distribution  of  Neutrons  from 
Be9 (a,n)B10  Reaction . --From  Shpetnyi.* 


*A.  I.  Shpetnyi,  Soviet  Physics  JETP,  5,  p.  357, 
October,  1957. 
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A  calibrated  Pu-Be  isotropic  neutron  source  with 
a  yield  of  3.5  x  106  neutrons  per  second  was  then  placed 
in  the  same  location  that  the  beryllium  target  assembly 
previously  occupied  and  the  neutrons  from  this  "known" 
source  were  counted  for  the  same  time  interval  using  the 
long  counter.   The  average  of  the  measurements  taken  at 
various  angles  using  the  Be9(d,n)B10  source  was  assumed 
to  be  the  measurement  due  to  an  "equivalent  isotropic" 
neutron  source.   This  average  value  was  then  divided  by 
the  measurement  taken  using  the  calibrated  source  and  the 
resultant  quotient  was  multiplied  by  the  strength  of  the 
calibrated  source  in  order  to  determine  the  Be9(d,n)B10 
source  strength  in  neutrons  per  second.   Thus,  the  "equiva- 
lent isotropic"  source  strength  for  various  deuteron  beam 
currents  was  determined.   The  results  are  shown  graphically 
in  Figure  III-3.   The  source  strength  for  a  beam  current 
of  1.5  microamperes  was  found  to  be  92.2  x  106  neutrons 
per  second . 

B.   Foils 


As  previously  mentioned  in  Chapter  II,  the  in 
herent  operating  characteristics  of  the  Van  de  Graaff 
accelerator  used  to  create  the  external  neutron  source 
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dictated  the  use  of  dysprosium  for  the  neutron  detection 
material.   The  stable  Dy164  isotope  useful  in  this  experi- 
ment comprises  28.2  per  cent  of  natural  dysprosium.   The 
other  isotopes  contribute  essentially  nothing  to  the  activ- 
ity measured  during  this  project. 

Since  no  information  concerning  the  self -absorption 
characteristics  of  dysprosium  could  be  found  in  the  litera- 
ture, the  variation  of  beta  counting  rate  with  foil  thick- 
ness was  determined  experimentally  following  standard 
procedures.11   Dysprosium  foils  of  several  different  thick- 
nesses were  irradiated  in  the  subcritical  reactor  at  the 
same  location  for  identical  time  intervals.   After  irradi- 
ation of  each  foil,  15  minutes  were  allowed  to  elapse  prior 
to  counting  so  that  the  activity  due  to  Dy165m  would  vanish. 
Each  foil  was  then  counted  for  a  period  of  10  minutes.   The 
results  are  shown  graphically  in  Figure  III-4.   This  self- 
absorption  curve  shows  a  broad  maximum  for  foil  thickness 
between  90  mg/cm2  and  130  mg/cm2.   This  corresponds  to 
thickness  between  0.0055  inches  and  0.0065  inches,  so  the 
foils  were  fabricated  with  thicknesses  as  close  as  possible 
to  6  thousandths  of  an  inch. 


X1C.    W.  Tittle,  Nucleonics,  8,  Ko •  6>  P«  5>    June, 
1951. 
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Thousands  of 
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Figure  III-4 


Variation  of  Beta  Counting  Rate  with  Foil 
Thickness  for  Dysprosium  Foils  1.93  cm  in 
Diameter . 
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The  foils  used  in  this  experiment  were  fabricated 
from  a  sheet  of  commercially  pure  natural  dysprosium  metal. 
Pieces  of  the  metal  were  turned  between  centers  on  a  lathe 
to  a  diameter  of  1  centimeter  and  were  then  ground  by  hand 
to  a  thickness  of  about  6  thousandths  of  an  inch. 

Upon  completion  of  the  cutting  and  grinding,  the 
foils  were  measured  as  accurately  as  possible.   The  foil 
weight  was  measured  on  a  Win.  Ainsworth  and  Sons  arm  balance 
to  an  accuracy  of  0.1  milligram.   The  diameters  were  deter- 
mined by  averaging  five  measurements  taken  with  a  micro- 
scope designed  for  measuring  tracks  made  by  particles  in 
nuclear  emulsions.   The  area  was  then  computed,  and  the 
foil  thickness  was  determined  in  units  of  milligrams  per 
square  centimeter.   The  results  of  these  measurements  and 
calculations  are  tabulated  in  Chapter  IV. 

C  .   Beta  Counting  Equipment 

A  liquid  scintillator,  composed  of  4  grams  per 
liter  p-terphenyl  and  0.1  gram  per  liter  "Popop"  dissolved 
in  toluene,  was  used  in  order  to  create  4ir  geometry  for 
beta  detection.12   This  scintillator  was  contained  in  a 


12 


T.  T.  Doss,  Unpublished  Thesis,  University  of 


Texas,  June,  1961. 
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cylindrical  aluminum  sleeve  with  a  removable  aluminum  cover. 
This  sleeve  (see  Figure  III-5)  screwed  into  a  machined 
brass  sleeve  which  covered  the  6810-A  photomult iplier  tube 
used  to  detect  the  scintillations  created  by  passage  of 
beta  particles  through  the  liquid  scintillator. 

Two  teflon  "0"  rings  were  used  to  create  a  seal 
between  the  photomult iplier  tube  and  the  aluminum  sleeve, 
thus  holding  the  liquid  scintillator  in  the  upper  chamber. 
The  removable  cover,  which  provided  easy  access  to  the 
liquid  scintillator,  had  to  be  designed  so  that  the  system 
was  entirely  light  tight  when  closed.   This  was  accomplished 
by  means  of  a  concentric  series  of  "tongues  and  grooves" 
machined  in  the  surfaces  of  the  aluminum  sleeve  and  the 
cover.   The  foils  were  suspended  in  the  liquid  scintillator 
by  means  of  a  small  spring  clip  attached  to  the  removable 
cover. 

The  dimensions  of  the  aluminum  sleeve  were  such 
that  the  amount  of  liquid  scintillator  surrounding  the  foil 
was  sufficient  to  absorb  all  the  beta  particles  emitted  by 
Dy165  decay.   The  maximum  range  of  the  most  energetic 
(1.3  Mev)  beta  particle  emitted  by  Dy165  was  determined 
by  the  equation13 


1 3  S  .  Glasstone^  Sourcebook  on  Atomic  Energy  (2d 
ed.,  Princeton,  N.J.:   D.  Van  Nostrand  Co.,  1958}.,  p.  187, 
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Figure  III-5 


Details  of  4tt  Geometry  Liquid  Scin' 
tillator  Beta  Detector. 
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Rmax(g/cm2) 


1.9 


for  0.5  Mev  <C  Emax  <  3  Mev 


For  Emax  =  1.3  Mev,  this  equation  gave  a  maximum  range  of 
0.595  g/cm2.   Toluene  has  a  density  of  0.866  g/cm3,  so 
Rmax  was  0,688  cm  in  the  liquid  scintillator.   The  smallest 
thickness  of  scintillator  adjacent  to  the  foil  was  0,750 
cm,  so  all  beta  particles  were  absorbed  prior  to  reaching 
the  walls  of  the  container. 

The  6810A  photomult iplier  tube  was  mounted  in  a 
Texas  Nuclear  Corporation  tube  base  utilizing  a  417A  cath- 
ode follower.   A  block  diagram  of  the  entire  beta  detect- 
ing and  counting  apparatus  used  in  this  experiment  is  shown 
in  Figure  III-6. 

The  proper  discriminator  setting  for  the  Baird- 
Atomic  amplifier  was  determined  by  increasing  the  setting 
from  zero  to  the  point  where  the  observed  counting  rate 
approached  the  background  level.   At  this  discriminator 
setting,  the  photomult iplier  tube  electronic  "noise"  was 
eliminated.   Undoubtedly,  detection  of  many  of  the  very 
low  energy  beta  particles  was  also  eliminated.   However, 
this  loss  is  corrected  by  the  FeFsa(^BC  ( d)  factor  previously 
discussed  in  Chapter  II. 
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The  factor  FeFsa£sc(d)  was  determined  experimen- 
tally by  coincidence  counting  techniques.   The  various 
foils  were  irradiated  long  enough  to  reach  saturation  for 
the  2.32  hour  half-life  activity  of  Dy165.   Fifteen  minutes 
were  allowed  to  elapse  between  the  time  of  removal  of  a 
foil  from  the  irradiating  flux  and  the  time  of  counting. 
A  foil  was  placed  in  the  beta  counting  apparatus  and  then 
a  gamma  ray  counting  apparatus  utilizing  a  Nal  crystal 
gamma  detector  was  positioned  adjacent  to  the  aluminum 
sleeve  which  covered  the  plastic  scintillator  for  the 
beta  counter.   Gamma  rays  emitted  co incidentally  with  the 
beta  decay  of  Dy165  passed  through  the  aluminum  sleeve 
and  some  were  detected  by  the  Nal  crystal. 

The  coincidence  count,  CR   ;  was  measured  first. 
The  gamma  ray  count,  C  ,  and  the  beta  particle  counts,  Cft 
and  Co  ( 4  Jt  counter)  were  then  measured  for  exactly  the 
same  time  interval,  and  the  observed  C   and  Co  values 
were  projected  back  along  the  proper  decay  curve  on  semi- 
log graph  paper  in  order  to  normalize  the  C  ,  Co,  and  C 
counts  to  the  same  time  interval.   The  factor  FeFsa&sc(d) 
was  then  determined  by  dividing  Cg^y  X  Cg(4it  counter)  by 
Cy    X  Cq  as  discussed  in  Chapter  II.   The  values  of 
FeFsa£_sc(d)  thus  determined  are  tabulated  with  other  foil 
data  in  Chapter  IV. . 
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D .   Experimental  Procedure 

The  foils  were  exposed,  one  at  a  time,  at  various 
locations  in  the  radial  access  port  of  the  reactor  assembly 
at  intervals  of  approximately  1.2  centimeters.   The  irradi- 
ation time  for  each  foil  was  5  minutes.   After  waiting  15 
minutes  after  removal  of  the  foil  from  the  reactor,  the 
beta  activity  from  the  decay  of  Dy165  was  counted  for  a 
period  of  10  minutes.   During  the  15  minute  waiting  period 
a  10  minute  background  count  was  taken  and  this  value  was 
assumed  to  be  the  background  count  during  the  foil  count- 
ing period. 

This  procedure  was  followed  for  both  the  cadmium 
covered  foil  measurements  and  the  bare  foil  measurements. 
The  same  foil  was  used  for  both  measurements  at  a  given 
location.   During  the  bare  foil  measurement,  a  polyethylene 
spacer  disc  was  used  to  replace  the  lower  cadmium  cover  in 
order  to  position  the  foil  in  exactly  the  same  location 
insofar  as  was  feasible.   The  actual  difference  in  location 
was  never  more  than  0.33  millimeters,  and  since  CC{^  is  only 
about  1  per  cent  of  0+^,  the  error  due  to  any  difference 
in  location  was  assumed  to  be  completely  negligible.   The 
cadmium  cover  used  was  in  the  form  of  a  cup  with  a  disc 
which  fitted  completely  inside  the  lip.   The  cadmium  used 
in  this  cover  was  0.051  centimeters  thick. 
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The  beam  current  (and  hence  the  source  strength) 
for  each  irradiation  was  determined  by  taking  20  beam  cur- 
rent readings,  equally  spaced  in  time,  during  the  5  minute 
irradiation  period.   The  average  of  these  20  readings  was 
the  value  used  for  the  beam  current.   All  foil  activity 
measurements  were  normalized  to  the  same  beam  current  of 
1.5  microamperes.   This  amounted  to  normalization  to  an 
external  source  strength  of  92.2  X  10s  neutrons  per  second, 
since  the  source  strength  is  directly  proportional  to  the 
beam  current,  and  since  92 ,2  X  106  neutrons  per  second  is 
the  source  strength  corresponding  to  1.5  microamperes  of 
deuteron  beam  current  (see  Figure  III-3 )  . 

The  scintillation  counter  was  kept  completely 
full  of  liquid  scintillator  at  all  times,  and  the  liquid 
scintillator  was  changed  frequently  in  order  to  prevent 
contamination.   The  foils  were  kept  as  clean  as  possible, 
and  were  washed  in  used  scintillator  each  time  prior  to 
being  inserted  into  the  detector. 

Prior  to  commencing  measurements  each  day  the 
counting  system  was  checked  for  proper  operation  using  a 
Na22  source  of  constant  activity.   This  source  gave  a 
counting  rate  of  28090  counts  per  2  minutes.   The  various 
equipment  settings  used  were  as  follows: 
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High  Voltage   

Amplifier  Discriminator 

Amplifier  Gain   .... 
Scalar  Discriminator 


1900  volts 

080 
6400 

0.0 


Measurements  were  taken  first  with  no  control  rod 
in  the  reactor  and  then  with  the  cadmium  control  rod  in  the 
reactor,  and  lastly  with  the  boron  control  rod  in  the  re- 
actor. The  control  rods  were  placed  along  the  longitudinal 
axis  and  covered  the  entire  length  of  the  axial  access  port 
of  the  reactor.  The  effect  of  the  boron  rod  was  small  and 
was  not  reported  in  this  thesis. 


CHAPTER  IV 


DATA    AND    GRAPHS 
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Figure  IV-1.   Radial  Neutron  Density  Distribution  (Experimental) 
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Comparison  of  Theoretical  and  Experimental  Radial 
Neutron  Density  Distribution  Curves  without  Con- 
trol Rod. 
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CHAPTER 


V 


RESULTS 


The  radial  neutron  density  distribution  curves 
as  determined  by  experimental  methods  are  shown  in  Fig- 
ure IV-1  for  the  reactor  with  and  without  the  control  rod 
in  position  along  the  longitudinal  axis  of  the  assembly. 
The  neutron  density  values  shown  here  have  been  normalized 
to  a  constant  source  strength  of  10   neutrons  per  second 
from  the  Be(d,n)  reaction.   Comparison  of  the  two  curves 
shows  a  large  depression  of  the  neutron  density  in  the 
vicinity  of  the  control  rod  and  a  general  decrease  through- 
out the  reactor.   This  is  to  be  expected  since  the  large 
number  of  neutrons  absorbed  by  the  cadmium  rod  are  no 
longer  available  to  diffuse  through  the  reactor  or  to  par- 
ticipate in  fission  events  which  release  other  neutrons 
to  the  system. 

The  single  velocity  group  diffusion  theory  tech- 
nique was  used  to  predict  the  neutron  density  distribution 
for  each  case  mathematically.   The  results  of  this  treat- 
ment are  shown  graphically  in  Figure  IV-2  for  the  reactor 
with  no  control  rod  and  in  Figure  IV-3  for  the  reactor  with 
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the  control  rod.   The  agreement  "between  the  theoretical 

and  experimental  dis tribu tions  could  be  termed  no  better 
than  fair  in  either  case. 

The  large  "bump"  observed  experimentally  at  the 
boundary  between  the  core  and  the  reflector  is  not  ade- 
quately represented  by  the  theoretical  dis tribu tions ,  al- 
though a  slight  change  in  slope  is  indicated  at  that  posi- 
tion.  It  is  interesting  to  note,  however,  that  the  areas 
under  the  theoretical  curves  are  essentially  the  same  as 
the  areas  under  the  corresponding  experimental  curves. 
This  implies  that  the  theory  does  provide  for  the  proper 
total  number  of  neutrons  although  the  distribution  is 
somewhat  distorted. 

It  is  believed  that  the  agreement  between  the 
di stribu tions  predicted  by  diffusion  theory  and  those  ob- 
served experimentally  would  have  been  better  if  a  more 
accurate  value  had  been  available  for  the  buckling  of  the 
reactor  used  in  this  experiment.   Calculations  show  that, 
if  the  value   B    =  0.02cm     were  used  in  this  case,  the 
neutron  density  predicted  theoretically  for  the  core- 
reflector  interface  would  have  been  almost  exactly  the 
same  as  that  experimentally  observed. 

Analysis  of  the  data  obtained  for  various  posi- 
tions at  the  core-reflector  interface  indicate  that  surfaces 
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of  equal  neutron  density  would  have  the  form  of  an  ellip- 
soid of  revolution  with  a  ratio  of  major  and  minor  semiaxes 
of  approximately  1.2  to  1.   This  implies  that,  if  the 
change  in  slope  at  the  core-reflector  interface  were  neg- 
lected, theoretical  calculations  for  the  system  might  be 
made  with  considerable  success  using  spherical  geometry. 
Calculations  involving  spherical  coordinates  are,  of  course, 
much  more  simple  than  are  those  involving  cylindrical  co- 
ordinates . 
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REACTOR    CONSTANTS    AND    CALCULATIONS 
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Core  Dimensions* 


Fuel 

Hole 

Element 

Thickness 

Diameter 

Diameter 

Volume 

Number 

(inches ) 

(inches) 

(inches) 

(cubic  inches) 

1 

1.900 

9.913 

1.282 

144.18 

2 

2.050 

9.908 

1.286 

155.38 

3 

2.070 

9.906 

1.280 

156 . 84 

4 

2.070 

9.905 

1.280 

156.84 

5 

2.045 

9.912 

1.280 

155.17 

6 

2.057 

9.934 

1.282 

156.76 

7 

1.035 

9.931 

1.280 

78.84 

8 

NOT 

USED 

9 

0.507 

9.920 

1.283 

38.52 

10 

NOT 

USED 

Total  Volume 

=  1,042.54 

Expressed  in  Terms  of  Cubic  Centimeters,  Centimeters,  and  Grams. 
Total  Volume  =   17,090.36  cm3 
Average  Diameter  =  9.914  in.  =  25.182  cm 
Total  Thickness  =  13.734  in.  =  34.884  cm 
Total  Weight  of  Core  =  18,603.2  grams 
Core  Density  =  1.08852  grams /cm3 
Total  Weight  of  U02  in  Core  =  2,668.49  grams 
Total  Weight  of  U235  in  Core  =  465.04  grams 
Total  Weight  of  U  in  Core  =  2,351.54  grams 


*Data  furnished  by  manufacturer,  Lockheed  Aircraft  Corp. , 
Marietta,  Ga. 
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REACTOR    CALCULATIONS 


14 


The  transport  mean  free  path,  A.tr,  for  a  medium  is 
given  by  the  equation 


tr 


Hv         *total  "  Zs  <2/3A> 


where 


Ztr  -  Ztotal    Zs  3A 


Z-  L  total    ^s    3A.J 


for  the  i  elements  of  which  the  medium  is  composed.   Thus, 
for  the  reactor  core, 


(c)     i  1 

'tr   -  "TcT  - 

L+  1.2792  cm  x 

tr 


=  0.78174  cm 


By  the  same  method  X         for  the  reflector  is 

tr 


(R) 


l 


i 


tr     ~TR7 

L  1. 2613  cm  1 

tr 


0.79283  cm, 


14Meghreglian  and  Holmes,  Reactor  Analysis  (New 
York:   McGraw-Hill  Book  Co.,  I960),  Chap.  5. 
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The  diffusion  coefficient,  D,  for  a  medium  is  given 
by  the  equation 


X 


D  = 


tr 


where  X-^r  is  the  transport  mean  free  path.   For  the  core  D 
is 


D„  = 


(c) 

tr 


0.78174  cm  =  0.26058  cm 
3 


and  for  the  reflector.  D   is 


D   = 


(R) 

tr     0.79283  cm 


=  0.26428  cm 


is 


is 


The  extrapolated  radius  p0  at  which  p    becomes  zero 


P0  ~    Pmeasured  +  2Dr  =  28.32  cm 


The  extrapolated  height  z0  at  which  p1    becomes  zero 


zo  ~  zmeasured  +  ^R  -  33.55  cm, 
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The  values  of  p  and  z  at  the  interface  "between  the 
core  and  the  reflector  are 

Pj_  =  12.59  cm 


and 


z  .  =  17 . 44  cm, 

l 


The  diffusion  length,  L,     is  defined  "by  the  equa- 


tion L2  =  D/Z  .   Thus,  for  th 


e    core 


2  D  0.26058    cm 

'c         "J^         0.0752    crn-1 


=  3.465  cm' 


and  for  the  reflector 


2     D      0.26428  cm     _  ___    2 
R     tRJ    0.0275  cm  x 
a 


l/L- 


The  constant  k2  is  defined  by  the  equation  k2 
Thus,  for  the  reflector 

(R) 


.2     I 

kR  ~  L 


R    DR 


—  =0. 1041  cm-2 


In    Chapter    I    it    was    determined    that    k       =    it/2z     .        Thus,    for 

r  z     '   o        ' 

the  reactor  used  in  this  experiment 
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k^    =  sA-  =  0.0468  cm  x 


Since  k2  is  defined  to  be  k2  =  k2  -  k2  ,  the  term 
R  R     p     z> 

k       can  be  determined  as  follows: 


k   =  Vk?  +  k2  =  -/O. 01084  +  0.00219  =  0.11415  cm  x 
p        ii  Z 


In  order  to  determine  kz,  the  equation 


it   PR        jt 
kz  tan  kzzi  -  2z   d   tan  gz   Zi 


o   c 


o 


must  be  solved  by  numerical  methods.   This  solution  yields 
a  value  of 


kz  =  0.0471  cm 


-l 


Using  the  value  of  B2  furnished  by  the  manufacturer 
of  the  reactor,  k  ,     can  now  be  determined  by  the  equation 


kp  =  Vb?  "  k? 


Since  B2  =  0.027172  cm  2,     this  relationship  gives  the  value 


of 
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k   =  0.1580  cm"1. 
P 


Evaluation  of  constants  involving  Bessel  Functions 


at  p   and  b   are  as  follows 


Kq(k  dJ    K.  3.233 

°  ,     P  °s  =  -2-, r  =  0.0045 

I  (k  p  )    I  (3.233 

o   p  o      o 


bQ  =  b  -  €  =  0.95  cm  -  0.71  cm  =  0.24  cm. 


Jq(Vo)  _  Jq(0-0379) 
Yo<Vo>    Yo(0.037y) 


=  -  0.465. 


For  the  case  with  no  control  rod,  where  A   = 
78.7  n/cm3  when  p  =  z  =  0,  the  constant  CR  is  evaluated  as 
follows : 


•B 


A, 


J0(kDPi) 


Ko(*pP0) 


=  7i 


Vv>i>  -  i0upPo)  Viyi) 

7  JQ(l.989) 

[K0(l.437l)  -  0.0045  I0(l.437l)] 


n/ 


cm 


iving  Cp  =  80.56  n/cm' 


When  the  control  rod  is  in  position,  A  is  taken 
to  be  the  maximum  experimental  value  in  order  to  correlate 
experiment  and  theory.   This  gives  the  value  Ac  =  54.3  n/cm3 
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Using  this  value,  CR  can  be  evaluated  by  the  relationship 


Jo(kpPi)  " 


CR  =  ACr" 


J  (k  b  ) 


o  \  p  o 


Ko(Vi)  - 


K°UPPo)  I  (.V) 
Io<V>o>   °   P  Pl) 


=  54.3 


[J  (1.989)  +  0.465  Y  (l. 989)1 
o  o        J 

[K(l.437)  -  0.0045  I  (1.437)] 


n/ 


cm 


giving  CR  =  112.73  n/cm3. 

In  summary,  the  equations  for  n(p,z)  with  all  the 
constants  evaluated  are  as  follows: 

(a)  Without  control  rod: 

n  (p,z)  =  78.7  J  (0.158p)  cos  0.0471z 

n  (p,z)  =  80.56[K  (0.1142p)  -  0 . 00451  ( 0 . 1142p )  ] 


R 


x  cos  0.0468z 


(b)  With  control  rod: 

n  (p,z)  =  54.3[J  (0.158p)  +  0.465Y  (0.158p)]cos  0.471z 

nR(p,z)  =  112.73[K  (0.1142p)  -  0 . 0045  I  ( 0 . 1142p ) ] 

X  cos  0 . 0468z 


These  equations  were  used  to  determine  the  theoretical  neutron 
density  distribution  curves  shown  in  Figures  IV-2  and  IV-3. 
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The  probable  error  of  a  function  F(x,y ,  z)  is  given 


fcy 


1  5 


p    MV 

x  axi 


p    MY2  +  fp    MY 
y  by)       V  z  ^z/ 


where  P.P.  and  P   are  the  probable  errors  in  x,  y,  and  z, 

x'  y'  z  '   '        ' 

respectively.   The  probable  errors  P.P.  and  P   are  defined 
j  *  x'      y'  z 

by  the  general  equation 


P   =  0. 6745cr 


vhere  u  is  the  standard  deviation  of  the  value  observed. 

The  function,  n(p,z),  with  which  this  experiment 
is  concerned,  is  related  to  the  measured  quantities  by  the 
equation 


n(0,»)  -  K  °th  fsp  w 


m  F   F   - 
e   saoi.sc 


(d) 


where  K  is  a  constant  having  the  value 


l5H.  Margenau  and  G.  M.  Murphy,  The  Mathematics 
of  Physics  and  Chemistry  (2d  ed. ;    Princeton,  N.J.:   D.  Van 
Nostrand  Co.,  1956),  Sections  13.32-36. 
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K 


1  x  10' 


2.2  x  105  x  38,412  X  1.011  X  92.2  x  106 


=  0.1269  x  10 


-5 


when  source  strength  of  1  x  10   neutrons  per  second  is  used 
The  probable  error  for  n(p,z)  can  thus  be  written 


Eq.  (1):   P   =  K  JfPn    - 

th  m  F„  F 


£. 


+  Pi 


'th 


e   sa&sc> 


sp  m  F   F 

e   sa&sc. 


°thFsP 


m   Fe  Fsa&sc> 


•(' 


2  ^i 

C   F      \  12 

PF  F       th  SP 

ersa&sc     2  V2 

*e    ^sa&scy 


assuming  K  has  no  error. 

The  data  obtained  experimentally  for  the  point 
p  =  12.62  cm,  z  =  0  will  be  used  to  determine  the  probable 
error  for  this  experiment. 

The  value  C^  is  obtained  by  the  relation 


Cth  =  Ctotal  "  Ccd 


where  C-tQ-ta±    and  CC(^  are  measured  quantities.   The  probable 

error  P_    is 
th 


'r    =  /  P2,       +  P2, 
uth       ^ total     ucd 
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where 


P2  =     [0.6745   \/l77475]2    =    81,134 


Ctotal 


and 


P2        =  [0.6745  y/!4:02]2    =  640 
Ccd 


Thus 


P         =  ^81134  +  640  =  286. 

cth 

The  probable  error  in  Fc,„(d)  is  called  P™    and 
use  of  the  average  error  of  1.5  per  cent  tabulated  by- 
Tittle16  for  the  function  F„-n(d)  gives  PF    =  0.01665. 

The  probable  error  in  the  foil  weight,  P  ,  is 

0.0001  grams. 

The  probable  error  Pp  p       for  the  self  absorp- 

e  sa&sc 

tion  and  self  scattering  correction  factor  is  determined 
by  the  relation 


16 
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F    F       o 
e    sa&s  c 


c^*  •  <%<«>  V 


-7       Cp°7 


'P      7y 


+  fpp     e  zj 


C.       C„(4n) 


C  CR(4jt) 

since  F  F   „    =  i  —2 — Z — 2 where  I    is  a  constant  havirif 

e  sa&sc        CpC7 

the  value  I    =    0.8962.   Using  the  values 


Pr       =    0.6745  a/c7 
ui 


in  the  above  relation  gives 


PF         =  0.000866 
' e  se&s  c 


Thus,  the  total  probable  error  is  determined  to 
be,  using  equation  (l), 

Pn ( experimental )  =  0.466  n/cm3, 


most  of  which  can  be  attributed  to  the  low  Co_7  count  rates 

experienced  when  experimental  determination  of  F  F  9       (d) 
*  *  e  sa&s  c 


was  accomplished,  and  to  the  probable  error  in  F 


sp 


This  value  represents  a  total  probable  error  of 

Pn ( experimental )  =  1.77  per  cent 

for  the  neutron  density  values  determined  in  this  experiment 

The  probable  errors  in  the  constant  terms  used  in 
the  development  of  the  expression  for  neutron  density  are 
those  give   by  the  authors  making  the  determinations  and 
are  : 


4-5 
r 
ny 


=  100  barns  for  a 


4-5 

ny 


=  800  barns 


Pa4-5m  - 
uny 


200  barns  for  a 


4-5m 

ny 


=  2000  barns 
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